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A practical real-time guidance algorithm has been developed for aerobraking vehicles that minimizes the post- 
aeropass A V requirement for orbit insertion while nearly minimizing the maximum heating rate and the mavimiim 
structural loads. The algorithm is general in the sense that a m inim um of assumptions is made, thus greatly 
reducing the number of parameters that must be determined prior to a given mission. An interesting feature is 
that in-plane guidance performance is tuned by adjusting one mission-dependent parameter, the bank margin; 
similarly, the out-of-plane guidance performance is tuned by adjusting a plane controller time constant. Other 
features of the algorithm are simplicity, efficiency, and case of use. The algorithm is designed for, but not necessary 
restricted to, a trimmed vehicle with bank angle modulation as the method of trajectory control. Performance 
of this guidance algorithm during flight in Earth’s atmosphere is examined by its use in an aerobraking testbed 
program. The performance inquiry extends to a wide range of entry speeds covering a number of potential mission 
applications. Favorable results have been obtained with a minimum of development effort, and directions for 
improvement of performance are indicated. 


Introduction 

I N previous years proposed aerobraking guidance concepts 1 were 
developed with computational efficiency as one of the major 
design goals. This was a carryover from a period when comput- 
ing capability was a fraction of what it is today. To accommodate 
this limited capability, simplifying assumptions were made in or- 
der to avoid numerical integration of the equations of motion. As 
a result, the guidance algorithms were not very general or adapt- 
able for future aerobraking mission planning, and laborious pre- 
mission analyses were required to determine appropriate values for 
mission-dependent parameters needed because of the simplifying 
assumptions. The speed of today’s computers permits more general 
real-time guidance modeling. 

Numerous aerobraking guidance strategies have been developed 
over the years. Fuhry 2 describes a guidance algorithm based on 
the numerical predictor-corrector of Gamble et al. 1 Fuhry’s 2 exit 
phase logic for apoapsis control is similar to the present in that 
the guidance algorithm determines the constant bank angle trajec- 
tory that takes the current vehicle state to the desired exit condi- 
tions. Shipley and Ward 3 discuss control strategies for aerobraking 
at Mars, providing a good discussion of issues relating to the tran- 
sition between entry and exit phases of flight; the present algorithm 
addresses this transition in a novel way. Willcockson 4 discusses the 
use of a numerical predictor- corrector algorithm that uses bank rate 
as a control variable for a continuously rolling vehicle, an approach 
not considered here. Braun et al. 5 discuss the use of a numerical 
predictor-corrector that uses a load relief strategy to explicitly limit 
heating rate during the entry phase. The present algorithm implicitly 
controls heating rate by issuing full lift-up commands (0 deg or 1 80 
deg bank angle commands, as appropriate for the vehicle) until the 
guidance logic determines that continuing to fly lift up may result 
in a premature skip-out. 

The algorithm described in this paper has been used in an aer- 
obraking simulation testbed program that allows the introduction 
of variations from nominal values of atmospheric density, vehicle 
mass and aerodynamic parameters, and on-board navigation hard- 
ware outputs. The initial conditions and target orbit can also be spec- 
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ified as desired. The results of many simulations conducted with this 
program have been reduced to a set of entry corridors depicted on 
the entry speed-entry flight-path angle plane. These corridors in- 
dicate the capabilities of the algorithm and suggest ways in which 
the two tunable parameters should be adjusted to obtain desired 
results. 

The paper is organized as follows. First an apoapsis controller is 
described. The strategy is to use constant bank angle trajectories to 
attain the desired apoapsis. Next a strategy is discussed that results in 
minimum A V and near-minimum peak heating rates and structural 
loading. The plane controller is then discussed. The orbit plane 
upon exit (inclination and longitude of ascending node) is implicitly 
controlled via the out-of-plane velocity. Finally, numerical results 
are presented in terms of entry corridors. 

Basic Apoapsis Controller 

The guidance algorithm examined in this paper is named the Gen- 
eralized Reusable Guidance for Optimal Aerobraking (GREGOAR) 
scheme. In this algorithm, a Newton-Raphson approach is used to 
determine the constant bank angle trajectory that will attain the 
desired apoapsis. Although the equations of motion are numer- 
ically integrated, for simplicity they are reduced to the follow- 
ing set of planar equations of motion 6 for a nonrotating, spherical 
planet: 

r = v sin y 

v = -D a - gsiny (1) 

y = L a cos<p/v + ( v/r — g/v)cosy 

Here r is vehicle position magnitude, u is inertial velocity magni- 
tude, y is the inertial flight-path angle, <j> is the bank angle, L a and 
D a are the aerodynamic lift and drag accelerations given by 

L a = [pv 2 r SC L lm D a = I pv?SC D /m (2) 
and g is the gravitational acceleration given by 

g = H/r 2 (3) 

In order to have adaptability with respect to variations in atmos- 
pheric density and with respect to uncertainties in the vehicle’s 
aerodynamic coefficients and mass, estimates of these quantities 
must be made during the aeropass. We assume that accelerometer 
measurements are available, so that measured lift and drag accel- 
eration magnitudes can be computed from the sensed acceleration. 
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Then, following Gamble et al., 1 the composite quantities pCi/m 
and pCofm can be estimated by 



(4) 


The estimates for each guidance cycle are input to a low-pass filter 
to obtain the values used in the guidance computations. 

The relative velocity v r in Eq. (2) can be computed by assuming 
the difference between inertial and relative velocity magnitudes is a 
constant throughout the trajectory. This is necessary because of the 
use of planar equations of motion. Assuming flight is hypersonic 
and that aerodynamic forces are significant over a small altitude 
range, the lift and drag coefficients can be modeled as constants, 
again for simplicity. 

The problem of determining the required bank angle <p is math- 
ematically that of determining the zero of the nonlinear constraint 
function 


= r at - r a (r, v, y)\ r = rcx (5) 

where r w is the radius at atmospheric exit, r at is the target radius of 
apoapsis, and r a is the actual radius of apoapsis at the simulated exit 
from the atmosphere. There are a number of techniques available 
for solving Eq. (5); the method used here is based on Newton’s 
method, with the Jacobian computed using a variational approach. 
To decrease the nonlinearity of the problem, the solution of Eq. 
(5) is obtained by iterating on the cosine of bank angle instead of 
on bank angle directly. If the solution for a given guidance cycle 
happens to result in |cos(<£)| > 1, which is referred to as guidance 
saturation, then the bank angle command is limited to 0 deg or 180 
deg as appropriate. For the next guidance cycle the limited solution 
from the previous guidance cycle is used as the initial guess. 

The set of variational equations used to solve Eq. (5) is developed 
here. The Newton-Raphson iteration from the cosine of required 
bank, c<f> y is written 

V 

c<p i+] = c<t>; - — —— (6) 

<M/d (c<t>) 


where 'I' is the constraint function (5) and d'i’/d (cip) is the Jacobian, 
both of which are evaluated at atmospheric exit Using conservation 
of energy and angular momentum, the constraint function can be 
expressed as 

W = 'l»(r, v, y) = 1 + 2/r(l/r 0 , - 1 fr)/v 2 - r 2 cos 2 y/r * (7) 

The Jacobian in Eq. (6) is expressed via the chain rule as 

dvp d'l' Br d y 9^ Bv d y 9^ d y 

d (c<t>) dr By d(c<t>) 9t> By d{c<j>) + By d (c</>) 


It is a routine matter to obtain analytical expressions for the partials 
of 'I' with respect to r, v, and y . The rest of the components needed in 
Eq. (8) must be obtained numerically as follows. Using the relation 

d (MQ/gg) = B(dx(t)/dt) 

dt da K } 


variational equations governing the motion can be written as 


d (Br/By) 
dr 


9 (dr /dr) 
By 


d(usiny) Bv 

= — siny -hucosy 

By By 


d (Bv/By) _ ^ (Bp Br/By t „9v/9y\ 

dt “ U °\JTr p ~) 

Br . 

+ 2— — sin y — g cos y 
r By 


d (By/B(c<f>)) = L a 
dt v 


(10) 


with the initial conditions 


zr~0 = to) = 0 
9y 

P-(t=k) = 0 ( 11 ) 


By 

B(c<!>) 


(t = 


to) = 0 


Equations (10) and (11) constitute an initial-value problem in or- 
dinary differential equations, which can be numerically integrated 
from the current time to atmospheric exit to give the required partials 
Br/By , Bv/By , and By/B(c<t>). In practice, the states r, v, y and the 
partials Br/By, Bv/By , By/B(c<j>) are propagated simultaneously to 
atmospheric exit, giving, in a single integration, all the information 
needed for use in Eq. (6). These equations are implemented in a 
noniterative mode; i.e., Eq. (6) is used only once per execution of 
the algorithm. The assumption used here is that the time between 
subsequent executions is small and therefore the vehicle state does 
not change significantly between executions. This eliminates the 
need for convergence criteria. 

An adaptive step-size fifth-order Runge-Kutta-Fehlberg numeri- 
cal integrator was chosen to perform the required integration because 
it is accurate and efficient and requires a relatively small amount of 
code and overhead. The integrated trajectories do not need to be 
computed very accurately because, due to navigation and modeling 
errors, the target can be attained with only limited accuracy anyway. 
Our experience in work on the aeroassist flight experiment (AFE) 
program has shown that 0. 1 km in position error in the integration 
is adequate. Typically, no more than 15-20 integration steps are 
required to compute a simulated trajectory. 

To ensure that the numerical integration process inside the guid- 
ance proceeds as intended, tests are performed after each integration 
step to detect simulated vehicles that are “trapped” in the atmosphere 
due to too much negative lift being modeled. If the (simulated) 
apoapsis is inside the sensible atmosphere or if the (simulated) min- 
imum altitude is below some specified value, say 50 km for Earth or 
30 km for Mars, then the numerical integration process is stopped. 
The amount of positive lift is then increased by incrementing the 
simulated cosine of bank by some specified amount, say ±0. 1 , and 
the integration process is started anew from t 0 - 

The equations and methodology described in the next section 
constitute the guidance logic used during the exit or ascent phase of 
the aeropass. This guidance phase is initiated at an appropriate point 
in the aeropass as determined by the entry phase logic discussed in 
the following section. 


Optimal Apoapsis Controller 

Extensive aerobraking trajectory optimization work has been 
done in recent years. 7-9 The strategy used here is suggested by 
the results of optimization work done by Miele and others, 10,11 who 
have shown two important results: 

1) If both atmospheric entry speed and flight-path angle are fixed, 
the optimal trajectory is a two-arc trajectory: an entry arc flown at 
full positive lift followed by an exit arc flown at full negative lift 

2) If only atmospheric entry speed is fixed, then the flight-path 
angle can be chosen so that the optimal trajectory is a one-arc tra- 
jectory flown at full negative lift. 

The trajectory of result 2 is a special case of the trajectory of 
result 1 where the switch to full negative lift occurs at entry. 

Here, “optimal” is used in the sense of minimum AV for 
postaeropass orbit insertion. 3,12 Optimization results have shown 
that the trajectories discussed here are characterized by near- 
minimum values of peak heating rate and structural loading. 12 
The optimal trajectories at aeropass exit satisfy the apoapsis con- 
straint (5). 

Numerical computation of this class of optimal trajectories is 
simple (once the form is known), in that, in the case of result 1, 
all that is needed is to iterate on a switching time (i.e., the time 
to switch from full positive lift to full negative lift) such that the 
vehicle attains the target apoapsis radius at exit. In the case of result 
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2, all that is needed is to determine the shallowest flight-path angle 
such that, at full negative lift, the vehicle attains the target apoapsis 
radius at exit. 

In a real-time guidance environment, result 2 does not have much 
practical use since the guidance discussed here becomes active at 
atmospheric entry and so obviously cannot influence an entry flight- 
path angle already determined by the exoatmospheric trajectory. 
Rather, result 2 is suggestive of a premission or pre-entry opti- 
mization (to determine entry flight-path angle) that must be done 
in conjunction with a real-time guidance algorithm based on result 
I. Thus the guidance strategy discussed here is derived with result 
1 as a basis. 

Result 1 cannot be used directly because the magnitude and ef- 
fects of dispersions (atmospheric, commanded bank angle rever- 
sals, navigation errors, etc.) that will be encountered during the exit 
(negative-lift) phase are unknown. If result 1 is used directly and, 
for example, the encountered atmospheric density were less than 
the atmospheric model predicted, there would not be enough lift 
capability to prevent a premature skip-out. What is needed is some 
performance reserve to account for the dispersions that the vehi- 
cle will encountered during the exit phase. This reserve is obtained 
by introducing a "bank margin," a premission-determined positive 
bank angle increment, which is used during the entry phase guid- 
ance mode in the simulation of the exit (negative-lift) arc. The exit 
phase arc is simulated as a constant bank trajectory flown at full 
negative lift less the bank margin (hereafter referred to as quasi-full 
negative lift). This has the effect of causing the guidance logic to 
start the exit phase earlier than it really needs to. Once the exit phase 
guidance mode has been triggered, it solves for the bank angle re- 
quired to hit the target apoapsis during each subsequent guidance 
cycle using the true vehicle lift capability (no bank margin included). 

Instead of solving for the switching time from entry phase to 
exit phase, the entry phase guidance computes one trajectory per 
guidance cycle simulating flight using quasi-full negative lift; in the 
meantime, it commands full lift-up. When the simulation indicates 
that a skip-out will occur despite being at quasi-full negative lift, 
the exit phase guidance is initiated. This strategy allows the vehi- 
cle to fly lift up as long as possible, thus minimizing penetration 
into the atmosphere and consequently reducing heating rates and 
structural loads while avoiding skip-out. It also maximizes perigee 
altitude, and so minimizes the AV needed to raise perigee out of 
the atmosphere after the aeropass is completed. 3 

As bank margin is increased, the guidance algorithm triggers the 
exit phase guidance earlier in the trajectory. This fact can be used in 
cases where the higher values of total heating loads associated with 
these optimal trajectories are more of a concern than peak heating 
rates. 

Note that all of the gains and many of the mission-dependent 
parameters so common in other guidance algorithms have been 
replaced in this optimal apoapsis controller by a single mission- 
dependent parameter, the bank margin. This feature of the algorithm 
makes it readily “portable" (or “reusable") from one vehicle to an- 
other and from one set of entry conditions to another, even to the 
extent of changing planets at which the aeropass is to take place. As 
a contrary example, in the scheme of Gamble et al. J gains on dy- 
namic pressure and altitude rate must be adjusted for every change 
in entry conditions, target conditions, and/or vehicle parameters. 

Another feature is that the time to transition from entry to exit 
phase is determined completely in real time: Prescribed trigger or 
transition velocities, which may have to be selected on the basis 
of trial and error, 1 3 are not necessary. A computationally desirable 
feature is that in the entry phase only one integrated trajectory need 
be generated (without the variational trajectories), and furthermore, 
early in the aeropass these are short-lived trajectories since the sim- 
ulated vehicles become trapped in the atmosphere. The exit phase, 
which is inherently more central processing unit (CPU) intensive, 
does not start until the time to exit has been reduced somewhat in 
the entry phase. 

Real-Time Orbital Plane Controller 

The form of the plane controller is based on a first-order sys- 
tem where the state variable being controlled is the out-of-plane 


velocity component, that is, the vehicle’s velocity component along 
a unit vector perpendicular to the desired orbit plane. Thus, if the 
magnitude of out-of-plane velocity is driven to a small value, the 
same is automatically done with the plane error. Starting with the 
homogeneous equation 


<*y y 

dt r 


= 0 


02 ) 


where y is the current out-of-plane velocity (measured positive along 
the unit vector antiparallel to the desired unit angular momentum 
vector) and r is a time constant, it is assumed that the out-of-plane 
lift component is the only out-of-plane force. Thus Eq. (12) may be 
written as 


L a sin0 + y/r = 0 (13) 

Solving for sin 0 and expanding L a , we obtain 


sin0 = (14) 

r PV?C L S } 

The sign of sin 0 gives the appropriate sign of the bank command 
whereas the magnitude of $in0 is a measure of the current plane 
error and thus can be used to determine when to switch the sign 
of the commanded bank angle. In practice, when the magnitude of 
sin <f> is larger than 1 , the sign of the commanded bank angle is set 
to that of sin 0: 


sgn (0 cmd ) = sgn (sin 0) (15) 

This commanded bank sign is used in subsequent guidance cycles 
until sin 0 changes sign again and |$in 0 1 becomes larger than 1. 

Inspection of Eq. (14) reveals that the magnitude of sin 0 becomes 
larger as dynamic pressure decreases. This results in increased sen- 
sitivity to plane error as atmospheric exit is approached and less 
sensitivity to plane error when dynamic pressure is large. This is 
precisely what is desired because larger plane errors deep in the 
atmosphere are acceptable; i.e., there is plenty of control authority 
available to correct for them; higher in the atmosphere less plane 
error is tolerated as exit is approached. 

Note that satisfactory guidance performance can be achieved by 
adjusting r. Thus, as in the apoapsis controller case, we have a 
situation wherein only one critical parameter value needs to be de- 
termined prior to a mission. An important point is that large values of 
r result in fewer roll reversals and large plane error at exit, whereas, 
on the other hand, small values of r result in more roll reversals 
and small plane error at exit. Thus, it is a matter of trading off roll 
reversals for plane error. 

Numerical Testing 

A version of the GREGOAR guidance logic discussed above has 
been coded into a FORTRAN subroutine. The code accepts a vari- 
ety of inputs, so that it may be used with different planets, vehicles, 
and entry conditions. A three-degree-of- freedom aerobraking guid- 
ance test bed program, the Aerobraking Orbital Transfer Technology 
(ABOTT) program, was used to examine the performance of GRE- 
GOAR. Inputs to ABOTT allow selection of vehicle mass, aerody- 
namic parameters, initial state vector, and other parameters neces- 
sary to conduct the simulation. In addition, multiple-run sets can 
be started with simulation initial conditions varied between runs, to 
allow statistical studies of the test algorithm to be performed. 

ABOTT assumes that the only control variable is the vehicle bank 
angle. Although the guidance models assume instantaneous bank an- 
gle responses to commands, ABOTT models finite bank rates and 
accelerations, with a maximum bank rate of 20 deg/s and a maxi- 
mum bank acceleration of 5 deg/s 2 (the same values used by Shipley 
and Ward 3 and in the AFE program 13 ). Angle-of-attack variations 
about a trim value specific to the vehicle are simulated by allowing 
the lift and drag coefficients to fluctuate in a Gaussian manner from 
one aeropass to another, with mean and standard deviation of the 
Gaussian specified by the user. Sideslip is assumed to be a constant 
zero. Variations in vehicle mass from one run to another can also 
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Table 1 Aeroassist flight experiment vehicle parameters 


Mass, 

kg 

Area, 

m 2 

Drag coefficient 

Lift coefficient 

Continuum Free molecular 

Continuum Free molecular 

1859.7 

14.314 

1.3549 2.06717 

-0.3807 -0.09574 


be applied. The mass is selected from a uniform distribution be- 
tween user-specified maximum and minimum values containing the 
vehicle nominal value. 

The program propagates the vehicle state from the initial time 
until either the vehicle exits the atmosphere or a preselected 
time limit is reached. At each time step within the propagation cy- 
cle the guidance algorithm is executed to compute the bank angle 
commands; auxiliary parameters such as atmospheric density, ac- 
celeration loads, aerodynamic heating, and dynamic pressure are 
calculated; and the equations of motion are numerically integrated 
using a fourth-order Runge-Kutta subroutine. When an aeropass 
simulation has been completed, summary data are output, including 
speed at entry, flight-path angle at entry, maximum aerodynamic 
load encountered, apogee altitude of the orbit at exit, and AV to 
circularize at the target altitude. 

For the multirun cases used in this study, ABOTT computes dis- 
persions to be applied to the vehicle mass, aerodynamic coefficients, 
and position and velocity vectors within the specified limits and se- 
lects an atmosphere model profile from a table specified by the user. 
Each aeropass simulation of a set is begun after the initial condi- 
tions have been modified in this way. When a sufficient number of 
runs have been made in the vicinity of a given initial point, i.e., a 
particular entry speed-entry flight-path angle combination, a step is 
made to the next initial point and a new set of runs is performed in 
its vicinity. Statistics for selected output parameters at each step are 
then compiled from the results. 

The vehicle parameters used in this study are those for the AFE 
vehicle. This vehicle's design parameters are well known and form 
a convenient data set 13 ; nominal values are given in Table 1. The 
stated lift and drag coefficients hold for the vehicle flying at its trim 
angle of attack of —17 deg. The combination of the AFE vehicle 
and GREGOAR guidance package will be referred to as the “study 
configuration/' 

In this study, atmospheric data are read from a table consist- 
ing of 101 cases generated by the Global Reference Atmosphere 
Model (GRAM) program. 14 The tabular data include density, pres- 
sure, temperature, and wind speed and direction as functions of 
altitude. The data from case 1 are nominal values and the remaining 
100 cases represent expected deviations about this nominal set. For 
a given run, ABOTT selects an atmospheric case index by using a 
random-number generator that chooses an integer from a uniform 
distribution between 1 and 101. 

Entry Corridor Definition 

The simulations are initiated at an altitude of 122 km above Earth. 
The entry speed is varied between 9.0 and 16.0 km/s in steps of 0.5 
km/s. At each step within this speed range, the entry flight-path 
angle is varied in steps of 0.01 deg between -3.0 deg and -10.0 
deg, so that the aerodynamic corridor 15 can be defined. The target 
condition is to exit the atmosphere on an orbit with an apogee of 
300 km in the same plane as the initial orbit. 

For the purposes of this study, an “entry corridor" is defined as 
that region on the entry speed-entry flight-path angle plane within 
which the study configuration can achieve the target orbit subject 
to selected constraints. For example, if the vehicle can perform a 
circularization bum of 150 m/s or less in AV, then a contour can be 
drawn in the above plane within which the A V to circularize at 300 
km would be less than or equal to this value. Regions in the plane 
outside this contour would violate the A V constraint and so would 
lie outside the entry corridor. Once the nominal entry speed for a 
mission is known and a nominal A V is defined for the vehicle, an 
entry corridor defined in this way would dictate the allowable range 
of entry flight-path angles. 

The upper boundary of an entry corridor is defined by the shal- 
lowest flight-path angles at atmospheric entry for which successful 


aeropasses can be performed. Entries at shallower angles will result 
in skip-outs. An upper boundary tends to be rather “hard," especially 
at the higher entry speeds, in that changes in entry flight-path angle 
as small as 0.01 deg can mean the difference between achieving the 
desired orbit and remaining on an escape trajectory upon exit from 
the atmosphere. The lower boundary of the corridor is defined by 
the sleepest flight-path angles for which the target condition can 
be reached. Steeper entries will result in suborbital trajectories or 
excessively low orbits. A lower boundary tends to be much “softer" 
than an upper one, since small changes in entry flight-path angle 
result in exit orbits with similar parameters. 

Corridors When All Variations Active 

All modeled sources of variation and error were assumed to 
be active when performing the runs described below. For each 
run the atmosphere profile was randomly selected, as described 
earlier. Since it was assumed that the vehicle mass would be 
fairly well known, it was allowed to vary within a range of ±1% 
from the nominal mass. The lift and drag coefficients varied in 
Gaussian fashion with distribution means centered at the nomi- 
nal AFE values and standard deviations of 0.1. The IMU error 
structure employed also reflected values characteristic of the AFE 
vehicle. 

Figure 1 plots the mean plus 2-sigma values of A V to circularize 
at 300 km vs the entry flight-path angle for an entry speed of 12.0 
km/s. For this cross section and for the complete entry corridor to 
be presented below, the GREGOAR bank margin was selected to be 
45 deg. Fitted curves like that in Fig. 1 were used to determine 2- 
sigma circularization A V contours on the entry corridor map. If one 
desired to specify a “nominal" circularization A V for an entry speed 
of 12.0 km/s, such a value could be found from Fig. 1 by selecting 
an acceptable entry flight-path angle and reading the corresponding 
A V off the curve. 

Figure 2 presents the entry corridor map obtained using the fitted- 
curve method mentioned above. Circularization AV contours at 
125 m/s (darkest shading), 150 m/s (medium shading), and 175 m/s 
(lightest shading) are plotted within a wedge-shaped boundary that 
represents the study configuration capability unconstrained by ve- 
hicle limitations on heating rate, structural load, etc. These contours 
were determined from many thousands of cases where all sources 
of variation were allowed to vary. 

Considering the 125-m/s zone, if a minimum 2-sigma corridor 
width of 0.5 deg is required to allow for navigation uncertainties, 
the highest speed for which the study configuration can be used is 
approximately 9.8 km/s. If this corridor width requirement is en- 
forced, the study configuration is marginally suited to transfer mis- 
sions from geosynchronous to low Earth orbit but is inadequate for 
higher energy missions. Considering the navigation infrastructure 
available in near-Earth space, a corridor width requirement of 0.5 
deg is probably unnecessarily stringent. For entry at Mars such a 
requirement may be necessary. 

If the A V budget is raised to 150 m/s, the entry speed that can 
be accommodated rises to at least 1 1.5 km/s. A bifurcation of the 
1 50-m/s corridor occurs at entry speeds above this point. The upper 
branch is at first about 0.2 deg in width, and the lower one is about 
0.15 deg in width. Thus neither branch meets the 0.5 deg minimum 
width requirement. 

If a 2-sigma A V budget of 1 75 m/s is allowed, the study configura- 
tion can accommodate entry speeds up to approximately 13.7 km/s. 
The corridor bifurcates at slightly higher speeds. Shortly after the 
split, neither branch meets the minimum-width constraint, though 
both are wider than the branches in the previous case. 

A further constraint that would quite likely be imposed on manned 
missions would be a maximum deceleration load limit. Previous 
studies 16 have proposed a value of 5 g for such a limit. Application 
of this constraint imposes the most severe limits on the mission 
suitability of the study configuration, as can be seen in Fig. 3. If 
the 2-sigma AV limit is 175 m/s, the 0.5 deg width requirement 
will make entry speeds above 12.5 km/s unacceptable due to the 
deceleration loads encountered. Even if this width requirement is 
dropped, the 5 -g constraint cuts off the 2-sigma corridor entirely by 
about 13.4 km/s. 
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Entry Flight Path Angle, Degrees 


Fig. 1 Entry corridor cross section at 12.0 km/s. 



Fig. 2 Entry corridor for all variations active. 


A similar corridor cutoff is seen in Braun et al. 5 at 13.9 km/s. but 
for a value of L/D of 0.1. For their case with an LfD of 0.3. which 
is similar to ours, the ^-limited corridor extends well past 14.0 km/s. 
A point to consider in this comparison is that the corridors presented 
in our paper result from variational cases, with the boundaries deter- 
mined by points 2-sigma from the mean of the result; if only nominal 
or mean values of atmospheric density and vehicle parameters had 
been used, the corridors would have been wider. This seems to have 
been the case for the corridor determinations of Braun et al., 5 which 
suggests one reason for their better corridor width results. When 
they performed an atmospheric dispersion analysis, they selected a 
nominal entry speed and flight-path angle combination (equivalent 
to selecting a “trough” point in Fig. 1) that resulted in a nominal 
value of circularization AV and then studied variations about this 
value. Their nominal A V value of 124.2 m/s at an entry speed of 14 
km/s compares favorably with our 2-sigma values at 12 km/s seen 
in Fig. 1, which already incorporate dispersion effects. 

We must admit that the comparative simplicity of our algorithm 
may also have contributed to the narrowness of its corridor at higher 
speeds. A technique to mitigate the narrowing of the corridor due to 


an acceleration limit involves using bank angle modulation to pro- 
vide load relief and so extend the flyable corridor to higher speeds. 15 
GREGOAR does not incorporate this strategy. 

Effect of Bank Margin on Entry Corridor 

One of the adjustable parameters of the GREGOAR guidance 
package is the bank margin. To examine the way in which changes 
in this parameter affect the entry corridor, two values besides the 
standard 45 deg were tried, 30 deg and 60 deg. A 2-sigma corridor 
cross section at 12.0 km/s for all three bank margin numbers is 
presented in Fig. 4. 

Recall that the bank margin is used to reduce the maximum lift 
the algorithm considers to be at its disposal when propagating the 
simulated trajectories prior to transition to the exit phase of flight. 
The larger the margin, the less lift considered, so that the transi- 
tion will be commanded earlier. Consequently, penetration into the 
atmosphere will be shallower at the transition time, and speed at 
transition time will be greater. This results in greater propulsive 
braking requirements to achieve the target orbit. 
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Entry Flight Path Angle, Degrees 

Fig. 4 Effect of bank margin on entry corridor cross section. 


As can be seen in Fig. 4, the curve for 60 deg bank margin is time constant, and examine the results. Figure 5 shows a corridor 

shifted to higher circularization A V values when compared with cross section at 12.0 km/s for the two time constant settings, where 

that for 45 deg margin. Acceptable entry flight-path angles are also the corridor is expressed in terms of the wedge angle upon exit. The 

shallower for the 60 deg margin setting. The corridor for this setting wedge angle combines the errors in inclination and right ascension 

is somewhat narrower and barely avoids bifurcation for an allowed of the node of the actual orbit with respect to the target orbit. In these 

AV of 175 m/s. If only 150 m/s is in the propulsion budget, the simulations the target orbit plane was taken to be identical to the 

corridor is closed for this margin setting. initial orbit plane, as stated earlier. In this figure the lower bounds of 

Considering a margin setting of 30 deg, fuller use of the vehicle’s each zone represent mean values of the exit wedge angle resulting 

lift will be made by the algorithm, transition will be delayed, and from their respective plane controller time constant values, whereas 

more of the vehicle’s velocity will be dissipated in the atmosphere. the upper bounds represent the mean plus 2-sigma exit wedge angle 

This can be seen in Fig. 4, where the curve for 30 deg margin reaches values. 

steeper entry angles, whereas minimum A V values are essentially Clearly, reducing the time constant reduces the plane error upon 

the same as for the 45 deg setting. Note that the trendency to bifur- exit: The minimum wedge angle achieved for a time constant of 30 

cation of the corridor is greatly reduced for the 30 deg value. This s was 0.069 deg. whereas the minimum for 45 s was 0. 1 7 1 deg. The 

is important for extending the corridor to higher entry speeds. The spread in expected error is also reduced, as seen by the much smaller 

5-e deceleration constraint may be somewhat mitigated as well. mean to 2-sigma band for the 30-s curves. This result confirms the 

' ' expectations described at the end of the derivation section. As stated 

Time Constant Effect on Corridor there, increased exit plane accuracy is purchased by an increase in 

The effect of the other adjustable parameter in GREGOAR, the the number of bank reversals during the aeropass. Recall that the 

plane controller time constant, remains to be examined. Since we bank reversals are not accomplished instantaneously, but are subject 

desire to improve the study configuration's performance, we select to a maximum rate of 20 deg/s and a maximum acceleration o 

a bank margin setting of 30 deg, use values of 45 and 30 s for the 5 deg/s 2 . In selecting a time constant, therefore, consideration must 
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be given to attitude control system requirements in performing bank 
reversals as well as the desired accuracy of the exit plane. Typically, 
four or five roll reversals are required with a time constant of 45 s, 
whereas five or six roll reversals are required with a time constant of 
30 s. The number of required roll reversals can probably be reduced 
with slightly more sophisticated logic. 

Conclusions 

The GREGOAR algorithm has shown itself to be resilient to vari- 
ations in flight parameters over a wide range of entry speeds. The 
tendency of the current version of this scheme to form bifurcated 
entry corridors could well limit its use on high-energy missions if 
large values of bank margin prove to be necessary . Further develop- 
ment work on the algorithm may reduce this characteristic splitting. 
In its current form the algorithm can be easily adapted for use in 
any number of aerobraking vehicles. Its minimum number of inputs 
and speed of execution, in addition to its robustness, recommend it 
highly for future use in aerobraking applications. 
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